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O . Abstract 

o 

A unified approach, for solving a wide class of single and many-body quantum prob- 
lems, commonly encountered in literature is developed based on a recently proposed 
method for finding solutions of linear differential equations. Apart from dealing with 
exactly and quasi-exactly solvable problems, the present approach makes transparent 
various properties of the familiar orthogonal polynomials and also the construction of 
their respective ladder operators. We illustrate the procedure for finding the approxi- 
mate eigenvalues and eigenfunctions of non-exactly solvable problems. 
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1 Introduction 
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The familiar series solution approach JlJ is routinely employed, when one encounters 
differential equations in quantum mechanical eigenvalue problems. This approach, not 
only is tedious to implement, but also throws little light on the underlying symmetries 
CU °f the problem at hand. Hence, one looks for alternate simpler methods, a classic 

example being the elegant ladder operator approach to the harmonic oscillator problem. 
Another often used approach is the factorization technique pioneered by Schrodinger 
§, Infeld and Hull §. 

In a number of cases, the underlying symmetries of the equation under considera- 
tion, have led to algebraic approaches based on group theory [|| and supersymmetry 
||. Most of these methods fail to generalize to many-body quantal systems and are 
| also not easily applicable to non-exactly solvable problems. 

In this paper, we give a pedagogical description of a recently developed method 
for solving linear differential equations || [7) and apply it to a number of single and 
many-body quantum problems. The proposed method is simple and assumes no special 
symmetry of the differential equation (DE) under consideration. The symmetries and 
the algebraic structure of the solution space emerge in a natural manner Q. 

The paper is organized as follows. In the following section, we give a brief descrip- 
tion of the proposed method of solving linear DEs, both, for the single and the many 
variable cases. For the purpose of illustration, we then consider the familiar harmonic 
oscillator problem, as also the related Hermite DE. Novel expressions for the solutions 
of the confluent hypergeometric and hypergeometric equations are then provided for 
later use. 
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Section 3 is devoted to the construction of the ladder operators. It is shown that the 
novel form of the solutions, in conjunction with the Baker-Campbell-Hausdorff (BCH) 
formula, lead to a straightforward construction of the ladder operators, much akin to 
the harmonic oscillator case. We then proceed to the quasi-exactly solvable problem in 
section 4, wherein, the utility of the present approach for finding approximate eigen- 
values and eigenfunctions for non-exactly solvable problems is also illustrated. Section 
5 is devoted to the many-body correlated systems, which are currently under intense 
study. The procedure and the subtleties involved in dealing with such multi-variate 
Hamiltonians are explicitly pointed out. We then conclude in section 6 after pointing 
out other areas, where this approach may find fruitful application. 

For the purpose of facilitating comparison with the standard literature, units have 
been appropriately chosen in different sections. 



2 A simple approach to familiar differential 
equations 

For simplicity, we will first consider the case of single variable linear DEs and point 
out its multi-variate generalization later. A single variable linear DE, as will become 
clear from the examples of later sections, can be cast in the form 



[F(D) + P{x, d/dx)} y{x) = 



(1) 



i 00 a n D n and a n 's are some 



where, D = xd/dx is the Euler operator, F(D) = J2n 
parameters; P(x,d/dx) can be an arbitrary polynomial function of x,d/dx and other 
operators. The solution to Eq. ([!]) can be written as || [Tj, 



y(x) 




■11 



F{D) 



P(x, d/dx) 



(2) 



provided, F(D)x x = 0; here C\ is constant. Before proceeding further, we list two 
important properties of the Euler operator to be extensively used in the text. Euler 
operator is diagonal in the space of monomials i.e., Dx x = Xx x and other operators 
carry definite degrees, with respect to the Euler operator i.e., [-D,O d ] = dO d , where 
d is the degree of the operator O d . For example, [-D,x 2 ] = 2x 2 and [D, d 2 /dx 2 ] = 
—2d 2 /dx 2 . Using the above results, it is easy to see that, the operator 1/F(D) is well 
defined in the above expression and will not lead to any singularity, if P(x,d/dx) does 
not contain any degree zero operator. 

The proof of Eq. (|2|) is straightforward and follows by direct substitution [||. 
Alternatively, since F(D)x x = 0, equating F(D)x x modulo C\ and Eq. (|l|), one finds, 

C x F(D)x x . 



[F(D)+P{x,d/dx)]y(x 
Rearranging the above equation in the form 

F <°> [» + FlD) 
and cancelling F(D), we obtain 

1 



-P(x,d/dx) 



1 + 



F(D) 



P(x,d/dx) 



y{x) = C x F{D)x } 



y{x) = c x x } 



(3) 



(4) 



(5) 
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This yields 



y(x) = c x 



1 + F{D) P ( x i d / dx ) 
which can be cast in the desired series form: 



(6) 



y(x) 



Ca£( 

m=0 



F(D) 



P(x,d/dx) 



It is explicit that, the above procedure connects the solution y[x) to the space of the 
monomials x . This fact will be exploited in the later sections for obtaining the ladder 
operators and explicate various properties of the solution space. 

The generalization of this method to a wide class of many-variable problems is 
immediate. Using the fact that, F(D)X X = has solutions, in the space of monomial 
symmetric functions Q , where D = = J2i x i^i the solutions of those multi- 

variate DEs, which can be separated into the form given in Eq. (|l]), can be solved like 
the single variable case. As will be seen later, this procedure enables one to solve a 
number of correlated many-body problems. 

For illustration, we consider the harmonic oscillator problem. The Schrodinger 
eigenvalue equation (in the units, fi=LV=m=l) 



dx 2 



+ (2E n 



Yv, 



o 



can be written in the form given in Eq. ([!]), after multiplying it by x 2 

[D -l)D + x 2 (2E n -x 2 )]^ n = . 



(7) 



(8) 



Here, F(D) = (D - 1)D and the condition F(D)x x = yields, A = or 1. Using Eq. 
(0), the solution for A = is, 



Co EC- 1 ) 



^ m=0 



1 



Co 



1 



2! 



(D 
■x* + 



DD 



(x 2 (2E -x 2 )) 



[2£ ] 2 , (2! + [2£ ] 2 ) 4 (4! + (2!) 2 [2£ ] + (2!)[2£;o] 3 ; 



4! 



2!6! 



x v + 



(9.) 



Here ipo is an expansion in powers of x, whose coefficients are polynomials in Eq. The 
above series can be written in a closed, square integrable form, Co exp(— x 2 /2), only 
when = 1/2. Analogously, A = 1, yields the first excited state. To find the n th 
excited state, one has to differentiate the Schrodinger equation (n — 2) number of times 
and subsequently multiply it by x n to produce a F(D) = x n -£^ = Y\^Z^{D — I), and 
proceed in a manner similar to the ground state case. 

It is clear that, our procedure yields a series solution, where additional conditions 
like square integrability has to be imposed to obtain physical eigenfunctions and their 
corresponding eigenvalues. Once, the ground state has been identified and for those 
cases, where ifj(x) = i/jqP(x), where P(x) is a polynomial, one can effortlessly obtain 
the polynomial part, as will be shown below. Proceeding with the harmonic oscillator 
case and writing 



if> a (x) = exp( — —)H a (x) 
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one can easily show that H a satisfies 



1 d 2 

D ~ a ~W HJX) 



where a = E n — 1/2. The solution of the DE 



(10) 



H a (x) = C a ]T ("I) 

m=0 



1 d 2 



(11) 



(D - a) 2 dx 2 

yields a polynomial only when a is an integer, since the operator d 2 /dx 2 reduces 
the degree of x a by two, in each step. Setting a = n in Eq. (10), we obtain the 
Hermite DE and E n = (n + 1/2) as the energy eigenvalue. Below, we give the algebraic 
manipulations required to cast the series solution of Eq. ( [Til ) into a form, not very 
familiar in the literature. For the Hermite DE, F(D) = D—n and P(x, d/dx) = — | 
the condition F(D)x x = yields A = n, hence, 



H n (x) = C n £ (-1) 



m=0 



1 1 d 2 



(D — n) 2 dx 2 



Using, [D , (d 2 /dx 2 )} = —2(d 2 /dx 2 ) and making use of the fact 



ds e -< D - n) 



(D-n) 



(12) 



(13) 



we can write, 



1 1 



d 2 



1 1 



d 2 



2 (D-n) dx 2 
Hence in general, 



2 (D-n) dx 2 

1 1 d 2 

2 (D-n) dx 2 



1 d 2 



1 



2 dx 2 (D - n - 2) 

2 



2dx 2 



1 



(D-n- 4) (D-n- 2) 



(14) 



1 1 



d 2 



2 (D-n) dx 2 



--— ) TT 



2dx 2 J (—2i)" 
m! I 4 a!x 2 / 



(15) 



Substituting Eq. (|Tq) in Eq. (12), we obtain, 

oo 

ffn(x) = C n £(-l) m 
m=0 



m! \ 4 <ix 2 



(16) 



a result, not commonly found in the literature [10]. The arbitrary constant C n is chosen 
to be 2 n , so that the polynomials obtained can match with the standard definition [11]. 
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The algebraic manipulations shown above can be applied to the confluent hyperge- 
ometric DE 

d 2 N d 

+ (7- 



x 



dx 2 



dx 



a 



$(a;7;x) 







(17) 



to give 



®(a,j,x) 



■11 



r(7) 



r( 7 



a 



■ cxp 



dx 2 



dx 



(18) 



where the normalization has been chosen appropriately. Likewise, the solution to the 
hypergeometric DE 



d 



d 



d 2 



d 



dx 2 

can be written as 



x + (a + (3 + 1) x— + a/3 - x-j^ ~ 1— 



dx 



dx 



F(a,/3; 7 ;x) =0 



(19) 



F(a 8-rx) = (-l)-P V{a - mi) 



cxp 




. x 



(20) 



For convenience a table has been provided at the end, which lists a number of commonly 
encountered DEs and the novel exponential forms of their solutions. It is worth pointing 
out that, the solutions of confluent hypergeometric and hypergeometric DEs given 
above, are polynomial solutions, provided a and (3 are negative integers. It should be 
noticed that, unlike the conventional expressions, the monomials in the above solutions 
are arranged in decreasing powers of x. 

The series solutions for the same can be obtained by a simple modification of the 
DE. Multiplying Eq. (|l7]) by x, we get 



2 d 2 d 2 d 

dx z ax ax 



ax 



$(q; 7 ;x)=0 , (21) 



2_d_ 

dx 



ax. The requirement 



for which, F(D) = (x 2 -£j + ixj^j and P(x,d/dx) = 

F(D)x x = 0, generates two roots, A = 0, 1 — 7, yielding the two linearly independent 
solutions. For A = 0, we obtain the familiar confluent hypergeometric series 



a a(a + 1) x 2 



(22) 



Similarly multiplying the hypergeometric DE with x yields two roots, A = 0, 1 — 7 and 
A = solution gives rise to the well-known, Gauss hypergeometric series. Since a num- 
ber of quantum mechanical problems can be related to confluent and hypergeometric 
DEs [12], we hope that the novel expressions given above and in the table will find 
physical applications. 



3 Construction of Ladder Operators 

In this section, we first derive the ladder operators for the harmonic oscillator, making 
use of the exponential form of the solution for the Hermite polynomials derived earlier 
and then proceed to the Coulomb problem and the Laguerre polynomials associated 
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with it. Subsequently, we outline the steps required for generalizing these results for 
other cases. 

The advantage of the present approach, which connects the solution space to the 
space of the monomials, lies in the fact that, at the level of the monomials, the ladder 
operators can be constructed easily, which facilitates the construction of the same 
at the level of the polynomials and the wave functions. The raising and lowering 
operators for the Hermite polynomials are derived, starting respectively from x and 
d/dx at the level of the monomials. For example, xx n = x n+1 , when operated by 
exp[— (l/4)(d 2 /dx 2 )] = exp(— A), from the left on both the sides and after introducing 
an identity operator suitably, yields 



e- A 2xe A 2 n e- A x n 



2 n+l e~ A x n+l 



which leads to 



2x- 



dx 



H n (x) = H n+1 (x) 



In deriving the above relation, we have made use of the BCH formula, 
e- A Be A = B + [B,A] + ^[[B,A],A} + ... . 



(23) 
(24) 

(25) 



The construction of the lowering operator is analogous to the above procedure. 



From -^x r ' 



or 



nx 



n-l 



one obtains, by performing a similarity transformation, 



dx 



-A, 



-^-H n (x) = 2nH n _ 1 {x) 
dx 



(26) 



(27) 



The creation and annihilation operators for the harmonic oscillator problem follow 
after one more similarity transformation: 



e 2 [2x —\e 2 



dx 



1 

e 2 H n (x) 



or 



and 



x — 



2 n n\ 

_d_' 

dx 



I 2(n + l) 
2 n + l {n + 1)! 



e " 2 H n+1 (x) 



tp n (x) = ^2{n + l)ip n+l {x) ; 



d £l 
e 2 — e 2 
dx 



1 _a£ 

-e 2 HJx) 



2n 



2 n n\ 



leads to 



Here 



— +x 

dx 



2 n ' l (n - 1)! 
2ni(j n - 1 (x) 



e 2 iJ„_i(x) 



™ = \ZJ ex P(-y) F «( x ) • 



(28) 
(29) 
(30) 
(31) 
(32) 



Proceeding in an analogous manner, we derive the ladder operators for the Coulomb 
problem. The polynomial solution of the Laguerre DE 



d 2 . d 
x—t-k + [a — x + lj- — h n 



dx 2 



dx 



L a n {x) = 



(33) 



6 



can be written as, 



(-l) n 



nl 



d 2 . ^.d 

- x d^~ {a + 1) Tx 



. X 



(34) 



With B = [xd 2 /dx 2 + (a + l)d/dx] and proceeding as in the Hermite polynomial case 
we obtain 



x — 2x- 



d_ 

dx 



(n + l)L% +1 (x) . 



(35) 



To construct the lowering operator, one should start with a lowering operator different 
from d/dx at the level of the monomials, since a similarity transformation on this 
operator does not lead to a closed form expression. The simplest one is the operator 
B itself, since it commutes with e~ B . Acting B on the monomial x n , we obtain 



d 2 . d 
x— + (a + l) 



n(n + a)x 



n-l 



dx 2 v ' ' dx 

and performing a similarity transformation through e~ B , one gets 

L%(x) = -(n + a)L%_ 1 (x) 



d 2 . .d 

x d^ 2+{a + l) Tx 



(36) 



(37) 



To obtain the ladder operators for the radial wave functions, Rn(x) = x l e~? L^(x) 
of the Coulomb Hamiltonian, one proceeds in a manner similar to the oscillator 
problem. Explicitly, the action of the raising and lowering operator, respectively are, 



and 



x d_ 2 d_ /(/ + !) ' 

2 dx dx 2 dx x 



x d_ x d?_ 2 d 1(1 + 1) 
2 dx dx 2 dx x 



K(x)= -( n + l)K +1 (x) 



(38) 



(39) 



It is worth pointing out that the value of a = (21 + 1) can also be changed, since 

d ,„ 



Tx Ll{x) 



-l n /n!)exp 



^ d 2 _ , +2 ,d_ 
dx 2 dx 



dx 



x 



(40) 



(41) 



from where, one obtains the standard result, 

±L«(x) = -LZ±\(x) . 

It is clear that a composite operator can be obtained keeping n unchanged, while chang- 
ing the value of a. The forms of the ladder operators given above are not unique, one 
can construct more complicated operators starting from different ones at the monomial 
level. 

The exponential form of the solutions enables one to construct ladder operator 
for hyper geometric, confluent hypergeometric and other polynomials and functions 
[||. which can be used to construct similar operators for the quantum mechanical 
problems associated with these polynomials, e.g., Morse, Pdschl- Teller, Eckart and 
other potentials 12, 13]. 
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4 Quasi-Exactly Solvable Problems 



This section is devoted to the study of QES problems [14, 15 1. These problems are 
intermediate to exactly and non-exactly solvable quantum potentials, in the sense that, 
only a part of the spectrum can be determined analytically. These potentials have 
attracted considerable attention in recent times, because of their connection to various 
physical problems. 

We illustrate our procedure, through the sextic oscillator (in the units fi = 2m = 
u = 1), whose eigenvalue equation is given by : 



d 2 

Of 

--r-^ + ax + 72 



)ip(x) = Eip(x) 



Asymptotic analysis suggests a trial wave function of the form, 

ip(x) = exp(— bx 4 )tjj(x) , 

which leads to, 



-^ + 2^x^ + (a + ^)x 2 



ip{x) = Evp(x) 



(42) 



(43) 



(44) 



where x 6 term has been removed by the condition lQb 2 = 7. One notices that, the 
operator O = (a + 3^/t~)x 2 + 2^/jx 3 d/dx increases the degree of i/j(x) by two, if tp(x) 
is a polynomial. Confining ourselves to polynomial solutions and assuming that the 
highest power of the monomial in ip{x) is n, one obtains, 



a 



— = 2n + 3 



(45) 



after imposing the condition that O does not increase the degree of the polynomial. 
This is the well-known relationship between the coupling parameters of the quasi- 
exactly solvable sextic oscillator. Taking n = 4 and 7 = 1 for simplicity, and after 
multiplying the above equation with x 2 : 



D(D - 1) + Ex 2 + 8x 4 - 2x 5 



d 

dx 







we get, 



Mx) = Co E 



,m=0 



1 



D(D-l) 



x 2 E Q + 8x 4 - 2x l 



dx 



(46) 



(47) 



Modulo Cq, the above series can be expanded as, 



Mx) = l-E - + (E 2 - 16)- + (64£ - El)- + 



(48) 



The monomials having degree greater than four vanish provided, i?o = 0, ±8. It can be 
explicitly checked that, for these values of E , Eq. (El) is satisfied. The eigenfunctions 



8 



corresponding to these three values are given by, 



and 



exp(-— )[l + 4x 2 + 2x 4 ] 
exp(-^)[l - ^x A ] , 
exp(-^)[l-4x 2 + 2x 4 ] 



(49) 
(50) 
(51) 



This procedure generalizes to a wide class of QES problems [16]. 

Below, we demonstrate the method of finding approximate eigenvalues and eigen- 
functions for non-exactly solvable problems, using the well studied anharmonic oscil- 
lator as the example : 



d 2 

-r-K + ax 2 + /3x A 
dx z 



En 



lj) n {x) = 



(52) 



Proceeding as before, ipo(x) can be written as, 



M*) = c { E(-i) r 



> m=0 



1 



(D-1)D 



E x' 



ax 



Px ( 



. 1 



which can be expanded as, 



i>o{x) 



1 - § x2 + k\ ( 2a + E ») x '~h. ( 24/3 " {UaEo + e3 ° 



x 6 + 



(53) 



(54) 



Although a number of schemes can be devised for the purpose of approximation, we 
consider the simplest one of starting with a trial function ipo( x ) = exp(— /xx 2 
and matching it with i/jq(x). Comparison of the first three terms yields, 



vx 



Eo 
2! 



and 



EL 

2! 

[IV - 



EL 

3! 



2a E% 
4! + 4! ' 
(UE + El) 
30 6! 



(55) 



The resulting cubic equation in energy, El - E a = 3/3/2, leads to one real root and 
two complex roots. Choosing the real root on physical grounds, one obtains, 



En 



2 l/3 



a 



A 



+ 



A 



3 • 2V3 



(56) 



where A 



1/3 



The value of En, obtained in the 



40.5/3 + (1640.25/3 2 - 108a 3 ) 1 / 2 
weak coupling regime, matches reasonably well with the earlier obtained results [17]. 
An approximate ipo can be obtained from Eq. (|54"|). One can easily improve upon the 
above scheme by taking better trial wave functions. Similar analysis can be carried 
out for the excited states. The above expansion of the wave function may be better 
amenable for a numerical treatment. For example, an accurate numerically determined 
energy value can lead to a good approximate wave function. 
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5 Many-body interacting systems 



In this section, we will be dealing with correlated many-body systems, particularly 
of the Calogero- Sutherland (l8| and Sutherland type |1S[| . These models have found 
application in diverse branches of physics like fluid flow, random matrix theory, novel 
statistics, quantum Hall effect and others [20, 21]. A number of methods e.g., Lax pair 
[E2l, Bethe-ansatz techniques 



2J, |25|] have 



and SV-extended Heisenberg algebra 
been employed for studying these systems. 

We start with the relatively difficult Sutherland model, where the particles are 
confined to a circle of circumference L. The two-body problem treated explicitly below, 
straight-forwardly generalizes to N particles. The Schrodinger equation is given by (in 
the units fi = m = 1) 



1± * 

i=l 



7T 



L 2 



sin [ir(xi - x 2 )/L] 



(57) 



Taking, Zj = e 2mx j/ L and writing Tp x ({zi}) = Ui,i^j z '% 
equation becomes, 



■P/2 



Zjf 'J\({zi}), the above 



J2 A 2 + P Z - 1± ^(Di - D 2 ) + Eq — E\ 



i=i 



zi - z 2 



Jx{{zi}) = 



(58) 



where, A = Zi £-, E\ = 2(±) 2 E X , E = 2(±) 2 E and E = (f ) 2 /5 2 , is the ground- 
state energy. Here, J\({zi}) is the polynomial part, which in the multivariate case is 
the well known Jack polynomial ||. Here, A is the degree of the symmetric function 
and {A} refers to different partitions of A. J2iE 2 is a diagonal operator in the space 
spanned by the monomial symmetric functions, TO{a}, with eigenvalues J2i=i • A 
monomial symmetric function is a symmetrized combination of monomials of definite 
degree. For example for two particle case, there are two monomial symmetric functions 
having degree two. These are 1712,0 = x\ + x 2 and m^i = x\x 2 - Readers are referred 
to Ref. [9] for more details about various symmetric functions and their properties. 
Rewriting Eq. (|58|) in the form, 



£(A 2 - x?) + f-^riPx -d 2 )+e + ^\ 2 - E\ 



z\ - z 2 

one can immediately show that 



Mizi}) = 



(59) 



Jx({zi}) 



E(-!) n 

. 71=0 



1 



Ei(D 2 -Xf) V *i 



z 2 



D 2 ) + E + Y,^-E x 

i j 
xm\({zi}) . 



(6' 



For the sake of convenience, we define 

1 



and 



5 
Z 



P 



Z\ + z 2 
z\ - z 2 



~ \ 2 ) 

(D.-D^+Eo + Y^^-Ex 



(61) 
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The action of S on m\({zi}) yields singularities, unless one chooses the coefficient of 
m\ in Z m\({zi}) to be zero; this condition yields the eigenvalue equation 



Using the above, one can write down the two particle Jack polynomial as, 



■/a({*}) = £(-£)' 



n=0 



1 



;( 



Z\ + z 2 



(D X - D 2 ) - ^(3 - 2i)\ 



xmx{{xi}) 

Starting from 1712,0 = z i + z \i it is straightforward to check that 

S n m 2 , = -2{l3) n m\ fl for n > 1 . 
Substituting the above result in Eq. fl60|), apart from C2, one obtains 



4/3(zi +Z2) Z = 4/3m? j0 
1 



(4/W „) = -2/W 



(62) 



J 2 = m 2 , + E(-1) B (-2)C9) B Nlo 



\n=l 

(00 
EH 3 )' 
n=0 

2/3 2 

m2 -° + T^ mi -o 



m l,0 



(63) 



which is the desired result. The above approach can be easily generalized to the N- 
particle case. 

Another class of many-body problems, which can be solved by the present approach 
is the Calogero-Sutherland model (CSM) and its generalizations [^]. Proceeding along 
the line, as for the Sutherland model, one finds the eigenvalues and eigenfunctions for 
the CSM. 

The Schrddinger equation for the CSM in the previous units, is given by, 



d 



1 



■■' 1 N 2 9 2 N 

2§a^ 1 2^-' 7 ' • "2 § (* ,•,■)* 



(64) 



where the wave function is of the form, JT? 

if> n (x) = ifoP n ({xi}) = ZGP n ({xi}) . 



(65) 



Here Z = Ui<j{xi - %j f, G = exp (-^zZi^i } , 9 2 = P{P ~ 1) and P n {{xi}) is a 
polynomial. After removing the ground state, the polynomial P n ({xi}) satisfies, 



22 x ia- + E °- En - oE 
^ ai; 2 ^r* 



— -ffE 1 — 

dxi ^ , . (Xj 3Jj) 



P„({ii}) = , (66) 
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where E = \N + \(3N{N - 1). Defining 




i+3 



^ (x^ Xj) $Xi 



i a 



one can easily see, following the procedure adopted for the Hermite DE that 



P n ({x i }) = C n e- A ^m {n} ({x i }) , (67) 



where m^ n y({xi}) is a monomial symmetric function of degree n. The corresponding 
energy is given by 



E n = E + n 



(68) 



One can use the above procedure to solve many other interacting systems. 



6 Conclusions 



In conclusion, we have presented a novel scheme to treat exactly, quasi-exactly and 
non-exactly solvable problems, which also extends to a wide class of many-body inter- 
acting systems. The procedure was used for the construction of the ladder operators 
for various orthogonal polynomials and the quantum systems associated with them. 
The approximation scheme presented needs further refinement. It should be analyzed 
in conjunction with computational tools for finding its efficacy as compared to other 
methods. The many-body problems presented here have deep connection with diverse 
branches of physics and mathematics. The fact that the procedure employed for solving 
them, connects the solution space of the problem under study to the space of mono- 
mials, will make it useful for constructing ladder operators for the many-variable case. 
This will throw light on the structure of the Hilbert space of these correlated systems. 
Some of these questions are currently under study and will be reported elsewhere. 

Acknowledgements: T. S. thanks U.G.C (India) for financial support through 
the JRF scheme. 
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TABLE I. Some frequently encountered DE and their novel solutions. 

Below, the differential equations from top to bottom, respectively, are Hermite, Laguerre, 
Legendre, Gegenbauer, Chebyshev Type I, Chebysgev Type II, Bessel, Confluent Hypergeometric 

and Hypergeometric. 



Differential Equation 



F(D), D = x{d/dx) 



Solution 



d_ 

' dx 



rrA. — r> — 1 d , 
X dx n 2 IF 2 



'dx 



H n (x) = 

LZ = 



d 2 

' dx 2 



X 



2 £* + ^lx--^n+l)-£}P n (x)=0 



dx 2 

c 2 d2 

dx 2 



+ {2\ + l)x£ -n{2X + n) - 



d 2 

dx- 



:C*(x) = 



x2 £+ x £- n2 



d 2 

dx 2 



T n (x) = 



P 2 d 2 
dx 2 



3x£-n(n + 2) 



d 2 

dx 2 



Un(x) = 



J±u(x) = 

x(l - x )£ I + { 1 -[a + (3 + l]x)£ - af3 
xF(a, (3, 7, x) = 



(D — n) 

(D-n) 

(D + n + l)(D-n) 
(D + n + 2X)(D -n) 

{D + n){D-n) 
{D + n + 2){D-n) 

{D + u){D-u) 

(D + a) 

(D + a)(D + f3) 



H n (x) = C n exp 



.1 d 2 

idx 1 



.X" 



L%(x) = C n exp 



P n (x) = C„exp 
C*(x) = C„exp 



T n (x) = C„exp 
U n (x) = C n exp 



J±u(x) = C± u exp 
<J>(a,7,x) = C_Q,exp 



. d 2 

'dx^ 



(a + 1) 



_d_ 

dx 



.X 



1 d 2 ' 

2(D+n+l) ch? 

1 d 2 
2(D+rt+2A) ch? 



1 d 2 



2(D+n) da? 
1 d 2 
2(D+n+2) dx 7 



1 X 2 



.X 



„n 



:X' 
.X r ' 



2(D±u) 



.X 



X dx^ ' dx 



F(a,(3,j,x) = C_( a>/ j)X 



x exp 



[D+(a,/3)| + T 'dx 



.X 



-09,a) 



The solution to the DE [F(D) + P(x, d/dx)} y(x) = 0, is, 

y(x) = c x {E™=o(-i) m [pfaPfa "} xA ' P rovided F ( D ) xX = °- 
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